Shot-noise anomalies in nondegenerate elastic diffusive conductors 
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We present a theoretical investigation of shot-noise properties in nondegenerate elastic diffusive 
conductors. Both Monte Carlo simulations and analytical approaches are used. Two new phenomena 
are found: (i) the display of enhanced shot noise for given energy dependences of the scattering time, 
and (ii) the recovery of full shot noise for asymptotic high applied bias. The first phenomenon is 
associated with the onset of negative differential conductivity in energy space that drives the system 
towards a dynamical electrical instability in excellent agreement with analytical predictions. The 
enhancement is found to be strongly amplified when the dimensionality in momentum space is 
lowered from 3 to 2 dimensions. The second phenomenon is due to the suppression of the effects of 
long range Coulomb correlations that takes place when the transit time becomes the shortest time 
scale in the system, and is common to both elastic and inelastic nondegenerate diffusive conductors. 
These phenomena shed new light in the understanding of the anomalous behavior of shot noise in 
mesoscopic conductors, which is a signature of correlations among different current pulses. 



I. INTRODUCTION 



The trend of scaling down the spatial dimensions of 
conducting samples is a demanding issue for realizing 
advanced electron devices and at the same time a fasci- 
nating realm for investigating new physical phenomena. 
Mesoscopic systems are presently a subject of intensive 
research and, being at intermediate sizes between atomic 
and micrometric structures, are in the position to satisfy 
this issue. In particular, a lot of attention is devoted to 
mesoscopic diffusive conductors, which are characterized 
by a sample length, L, much smaller than the inelastic 
mean free path, l in , and much longer than the elastic 
mean free path, l e i. These systems exhibit shot noiseEJ 
when the applied voltage, V, is much higher than the 
thermal voltage Vr — ksT/q, where fc# is the Boltz- 
mann constant, T the bath temperature and q the unit 
charge. Shot noise can be conveniently studied in terms 
of the Fano factor 7 defined as 7 = Si (0) / (2ql) , where 
5/(0) is the low frequency current spectral density and 
/ the electric current. A value of 7 = 1 (full shot noise) 
implies the absence of any correlations among different 
current pulses, while 7 / 1 is a signature of the pres- 
ence of correlations. In this latter case, both suppressed 
( 7 < 1) and enhanced (7 > 1) shot noise are known to 
be possible. 

In degenerate mesoscopic conductors at T = Qua the 



correlations induced by the Pauli exclusion principle give 
rise to a universal Fano factor 7 = l/ipJ^^Qjotr^st, in 
nondegenerate mesoscopic conductor»»BBEl the 
Fano factor is in general a function of the applied bias 
and its value depends on the interplay between the long 
range Coulomb interaction and the energy dependence 
of the elastic scattering time. For the relevant case of a 
three dimensional conductor (in momentum space) with 
an energy dependence of the scattering timti-Qf the form 
r oc e Q , Monte Carlo (MC) simulationsoM have re- 
ported a systematic analysis of the Fano factor in the 
range of values —1.5 < a < 1. Results showed values 
of 7 close to 1/3 for a between and 1 and V rang- 
ing between 40 and 120 Vr- On the other hand, the 
values of 7 were found to increase systematically from 
1/3 towards 1 for values of a de^asies-from to —1.5. 
The existing analytical theorieaD'BmiJ'lLil provide a rea- 
sonable interpretation of these results in th e lim ited range 
of values -qL\A < a < 0.5, being^,bsentE2lEil or failing 
completelyQtlE -outside this range. tl In this context, two 
of the Authorai3 have recently pointed out the possibil- 
ity that the above results only provide a partial descrip- 
tion of the shot-noise properties of nondegenerate con- 
ductors. In particular, enhanced shot noise (i.e. 7 > 1) 
was predicted for given values of a, and full shot noise 
(with 7 = 1) was expected at asymptotically high bias for 
a > —3/2. We conclude that the scenario of shot noise 
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in nondegenerate diffusive conductors, while phenomeno- 
logical rich of interesting features, is still incomplete and 
waiting for a microscopic physical interpretation. 

The aim of this work is precisely to address this is- 
sue. To this purpose, we have performed extensive MC 
simulations for a nondegenerate elastic diffusive conduc- 
tor covering values of a, applied bias, and dimensionality 
that are necessary for a unifying and complete analysis 
of the physical problem, and that enable a valuable test 
of the theoretical predictionst3 to be carried out. The 
paper is organized as follows. Section || describes the 
physical model and the MC technique used in the cal- 
culations. Section [II summarizes the theoretical predic- 
tions. Section IV presents the results of MC simulations 



proving the possibility for nondegenerate elastic diffu- 
sive conductors to display enhanced as well as full shot 
noise. The comparison between elastic and nonelastic 
transport regimes is also carried out to outline an in- 
teresting anomalous crossover from thermal to full shot 
noise exhibited by diffusive conductors. Section [v] draws 
the main conclusions of this work. 



II. PHYSICAL MODEL 

Following previous works nondegenerate diffusive 
conductors are modeled as a simple structure consist- 
ing of a moderately-doped semiconductor active region 
of length L and constant cross-section A sandwiched 
between two heavily doped contacts that, by inject- 
ing carriers., into the active region, act as ideal thermal 



reservoirsJij According with the nondegenerate nature of 
the carrier injection and the thermal character of con- 
tacts, carriers are injected into the structure following a 
Maxwellian velocity distribution at the lattice tempera- 
ture T and a Poissonian time statistics, with injection 
rate Y = (1/2) n c Vth, where n c is the contact density and 
v t h = \j2ksT /itm the thermal velocity of injected carri- 
ers, with m being the carriers effective mass. Once carri- 
ers are injected, their dynamics (under the action of scat- 
tering) is simulated in the active region of the structure 
by means of an ensemble MC self-consistently coupled 
with a Poisson solver (with constant voltage conditions, 
at the boundaries) to account for Coulomb correlations 
If not stated otherwise, we assume that the scattering 
processes taking place in the active region of the sample 
are elastic and isotropic and that the energy dependence 
of the scattering time is of the form r(e) = r^e 01 . In prin- 
ciple, the exponent a can take any value, although only 
some particular values correspond to well known scat- 
tering mechanisms in a semiconductor model. As exam- 
ples, a = —1/2 corresponds to scattering with acoustic 
phonons by deformation potentials, a = to neutral im- 
purities, a = 1/2 to acoustic piezoelectric phonons-and 
a = 3/2 to short-range ionized impurity scatteringo All 
these scattering mechanisms are important at low tem- 
peratures, which represent the typical conditions for an 



experimental test of theoretical predictions. 

The structure is assumed to be sufficiently thick in 
transversal dimensions to allow for a one-dimensional 
electrostatic treatment. Accordingly, the Poisson equa- 
tion is only solved in the direction of the current flow. 
In computer simulations one can consider any value of 
the momentum-space dimensionality d. Thus, even if 
the physically relevant case is d = 3, we will analyze also 
the case d = 2, since for this value of d there exist an- 
alytical predictions based on the same assumptions on 
what concerns the caimli ng ;. o i j t he kinetic equation with 
the Poisson equation. TO'EMia As indicated in Ref. [ §, 
this case corresponds to a hypothetical "flatland" whose 
physical realization would be a layered material with each 
layer containing a two-dimensional electron gas. 

Due to their very high carrier concentration as com- 
pared to that of the active region, contacts are assumed 
to have no voltage drop inside them. The applied voltage 
between the contacts is then considered to be constant 
in time. As concerns noise calculations, this situation 
corresponds to current-mode operationJlj 

In our analysis of shot noise we are mostly interested in 
determining 5/(0), which, in MC simulations, is obtained 
from the time integration of the autocorrelation function 
of current fluctuations Ci(t). The calculation of C/(i) 
is performed from the time series I(t) provided by the 
MC simulation. J(t) is the instantaneous steady-state 
total current as measured in the outside circuit, which, 
for the assumed geometry, is also the same through each 
cross-sectional area of the device. For the semiclassical 
case considered here, the calculation of I{t) is performed 
starting from the definition of the microscopic single par- 
ticle distribution 



N{t) 

/(k,r,t) = ]T<5[k 



ki(t)]*[r -!■«(*)], 



(1) 



where N(t) is the instantaneous number of electrons in- 
side the sample, which implies the conduction current 
density 



j c (r, t) 



q J v/(k, r, t)dk = q ^(t)6{r - nit)} 



(2) 



where Vj(i) is the instantaneous group- velocity deter- 
mined through a band structure model in k-space. We 
notice that the above distribution function does contain 
fluctuations and has the following properties. j—Jts en- 
semble average satisfies the Boltzmann equatior£3 while, 
in the linear approximation of small deviations from 
its average stationary-value, it satisfies the Boltzmann- 
Langevin equation 



(-+S)/(k,r,t)=y(k,r,t), 



(3) 



with S the linearized scattering operator and i^(k, r, t) 
the fluctuating Langevin source related to SBa From 
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the solenoidal property of the total current-density, and 
noticing that the displacement current is implicitU^taken 
into account by constant-voltage conditions, it ia_3 



N(t) 



(4) 



1=1 



where Vi (t) is the instantaneous carrier velocity along the 
field direction, known from the MC simulation. We re- 
mark that under steady-state conditions I(t) is a stochas- 
tic quantity which enables the calculation of average val- 
ues and correlation functions of its fluctuations with re- 
spect to average values. The basic sources of fluctua- 
tions in our system come from the instantaneous changes 
of carrier velocity due to scattering mechanisms and 
of carrier number inside the device due to injection- 
extraction processes. Thus, velocity, number and their 
cross-correlation contributions to the total current fluc- 
tuations are naturally accounted for together with their 
coupling with the fluctuating self-consistent field. We 
remark also that the ensemble MC simulation, by ap- 
plying semiclassical dynamics to describe the motion of 
each particle between scattering events, is tantamount 
to a calculation of the microscopic distribution func- 
tion /(k, r, t) which, by a suitable discretization of phase 
space, could be collected during the simulation. Since 
two-particle interaction is here neglected (i.e. the corre- 
sponding Boltzmann equation is linear) and the collision 
rates are Markovian, the present ensemble MC simulation 
provides a numerical solution of the Boltzmann-Langevin 
equation (||). 

Here we are primarily interested in Ci{t). This quan- 
tity is directly calculated from the MC by the following 
procedure. After neglecting an initial transient, the val- 
ues of I{t) are recorded in a time grid of step size At along 
the total simulation time MAt, with M integer. Then, 
by defining the time length in which the correlation func- 
tion should be calculated as mAt, with to integer, Ci(t) 
is obtained from the algorithmli3 
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C I (JAt) = I(t')I(t' + jAt)-I 

1 M — 771 -. in 

T7 E J(iA*)7[(*+i)Ai] -hrrE 7 ^)] 2 ( 5 ) 



M 



i = l 



where the bar denotes time average (ergodicity is implic- 
itly assumed), and j = 0,1,..., to ; M > to. Typically, 
in our calculations M > 500 x m, to > 500 . The corre- 
sponding spectral density Si(f) is determined by Fourier 
transformation of the Cjit) so obtained. 

For the calculations we have used the following param- 
eters: T — 300 K , effective mass to — 0.25 mo, with too 
the free electron mass, e = 11.7 e (with e the vacuum 
permittivity), L — 200 nm and n c — 4 x 10 17 cto~ 3 . For 
these values L/Le, c = 30.9 ~S> 1. To test theoretical pre- 
dictions we have performed simulations for systems with 
d = 2, 3 and the following values of a (and associated 
values of t ) -2 (8.35 x 10~ 16 s eV 2 ), -3/2 (6.76 x lO" 16 



s eV 3 / 2 ), -1 (3.88 x 10~ 16 s eV), -1/2 (4.41 x lO" 16 s 
eV 1 / 2 ) and (2.00 x 10~ 15 s). The values of t are cho- 
sen in such a way to ensure diffusive transport through a 
sufficiently high number of scattering events (over 10 2 ) in 
the active region. Poisson equation is solved on a space 
mesh with 100 nodes. The number of simulated particles 
ranges between 100 and 2000 depending on a and bias 
conditions. Also depending on a, the time step used to 
solve Poisson equation ranges between 2 and 0.1 fs. The 
numerical uncertainty is confined within 10% in average 
and 20% at worst. 



III. THEORETICAL PREDICTIONS 

Before Before presenting the Monte Carlo results, 
we will briefly summarize and elaborate the theoretical 
predictions derived in Ref. [ [l^] concerning the shot 
noise properties of nondegenerate elastic diffusive con- 
ductors. In that work, by means of a formal treatment of 
the Boltzmann-Langevin-Poisson set of equations it was 
shown that the Fano factor for this class of systems can 
be exactly decomposed into the sum of three contribu- 
tions as: 



7 = Tin +14> + lin,<t> 



(0) 



The first contribution, ji n , is related to the intrinsic 
current fluctuations in the presence of a static non- 
sclfconsistent electric potential. The second contribution, 
70, is related to the current fluctuations induced by the 
fluctuations of the self-consistent potential. The third 
contribution, 7i W) 0, is related to the cross-correlations 
between the previous two sources of fluctuations. Un- 
der far from equilibrium conditions (i.e. V/Vr > 3), and 
for a nondegenerate conductor, the two following prop- 
erties are satisfied: ji n = 1, and 70 > 0. As a con- 
sequence of these properties, theory predicts the follow- 
ing three regimes of shot noise: (a) enhanced shot noise 
when 7^0 > -7^, (b) full shot noise when 70 = -7^,0, 
(c) suppressed shot noise when 7^0 < —70- Con- 
cerning enhanced shot noise, Eq.(j(|) and the properties 
ji n = 1, 70 > imply that 7.^0 > is a sufficient 
condition to obtain enhanced shot noise. This condition 
was found to be physically realized in the presence of 
a negative differential conductivity in energy space, i.e. 
cr'(s) — da(e)/dE n -f- j 1 where the conductivity in energy 
space is given by 



<j{e) = ^q 2 v 2 {e)v(e)T(e) 



(7) 



with v{e) = flm (2me/h 2 ^j d ^ 2 1 the density of states 
in a d— dimensional momentum space, e the kinetic en- 
ergy, to the effective mass, h the Planck constant, and 
n = 27r d / 2 /r(d/2) the surface area of the unit sphere in 
d dimensions with T(z) being the gamma function. More- 
over, v(e) = \j2ejm is the velocity and r(e) the elastic 
scattering time. Therefore, for r(e) = ro£ a it is 
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o-'(e) = C(a + d/2)e 



a+d/2-l 



(8) 



where C is a positive constant. We conclude that a < 
—d/2, by implying a'(e) < 0, represents a sufficient con- 
dition for the presence of enhanced shot »aise. This con- 
stitutes the main prediction of the theorytj which is here 
under test. 

Concerning full shot noise, Eq. (^) and the property 
7;„ = 1 imply that 7^ +7^^ = is a sufficient condition 
to obtain full shot noise since it corresponds to indepen- 
dent current pulses through the device. This condition 
can be realized by the two alternatives 7^ = 7i n ^ = or 
70 = ~lin,<j> ¥= 0. The first alternative, 74, = = 0, 
corresponds to washing out completely the effect of long 
range Coulomb correlations. It can be realized under 
three means: (i) on a space scale, through the condition 
that L < Ld c where Ld c — \J efcgT / q 2 n c is the contact 
Debye screening length calculated using the concentra- 
tion of carriers at the contacts, with e the static dielec- 
tric constant of the semiconductor and n c the contact 
concentration; (ii) on a time scale, through the condition 
t~t < Td, where tt is the transit time and T4 the dielectric 
relaxation time, when L > Lo a ', and (iii) on the energy 
dependence of the scattering time, through the condition 
for the energy exponent a = —d/2 (thus generalizing to 
any dimensionality the pioneer 3D result of NagaevEJ)..-, 

Condition (i) was tested by numerical simulationsu 
that showed full shot noise for V/Vr > 3 at Debye lengths 
comparable with the sample length. 

In condition (ii), being the sample length longer than 
the Debye screening length, there is room for Coulomb 
correlations to be effective. However, since the transit 
time depends on the bias, the onset of shot noise will start 
at voltages in general higher than those satisfying the 
condition V/Vr > 3. The bias value at which shot noise 
will appear can be roughly estimated from the equality 
tt = L 2 1 [iVshot = T~d = e/q/j,n c , with fi being an effective 
mobility. This gives V shot = qn c L 2 /e = (k B T / q)L 2 / L 2 D ^ 
which, since L > Lr> c , will lead to V s hot 3> Vr- There- 
fore, in this case we expect an anomalous transition be- 
tween thermal and shot noise for biases Vr < V < V s hot- 
It is worth noting that the expression for V s hot just de- 
rived coincides with the bias at which L = L s , where 
L s = y' eV/qn c is the length wed to characterize the 
space-charge limited conditions This makes the condi- 
tion for the appearance of full shot noise tt < Td equiv- 
alent to the concUtion—Ls < L, thus explaining why the 
existing theoriesJraBy which apply for L s > L , have 
failed in predicting the asymptotic presence of full shot 
noise. In any case, these theories are expected to remain 
valid for biases Vr < V < V s hot, and their predictions 
can be taken to explain at least partially the anomalous 
crossover between thermal and shot noise where 7 takes 
values below unity. 

In condition (iii), the crossover between thermal and 
full shot noise is monotonous and the Fano factor will 
not take values below one for any bias. 

As regards the second alternative to obtain full shot 



noise, 70 = -7^^ =^ 0, we have not found physically 
feasible conditions to accomplish it. From the previous 
discussion, it follows that there are two main situations 
which should be tested by a direct microscopic simula- 
tion, namely: the possibility to observe enhanced shot 
noise for a < —d/2 at different dimensionalities, and the 
possibility to observe full shot noise for L > Lr> and 
V > V shot . 



IV. MONTE CARLO RESULTS 

In this section we report the results of MC simulations 
with the objective of providing a complete picture of the 
shot-noise properties of nondegenerate elastic diffusive 
conductors and o f te sting the theoretical predictions for- 
mulated in Sec. 



III. Initially we will illustrate the dif- 



ferent conditions under which enhanced, suppressed, and 
full shot noise can be obtained in elastic diffusive con- 
ductors. Then, the crossover between thermal and shot 
noise will be detailed by comparing the noise power of an 
elastic diffusive conductor with that of an inelastic one. 



A. Enhanced, suppressed and full shot noise 

Figure |l| reports the Fano factor as a function of the 
normalized voltage for a major set of simulations per- 
formed with different values of the exponent a. Figure 
1(a) refers to the 3D case and Fig. 1(b) to the 2D case. 
The inset in each figure displays the I — V character- 
istics in terms of the saturation current I s = qTA = 
(l/2)qn c VthA, where A is the cross-sectional area of the 
structure. For voltages below 10 Vr the conductor ap- 
proaches thermal noise conditions, which explains a sys- 
tematic trend of 7 to values greater than unity. By con- 
trast, in the relevant region of voltages from 10 to 100 
Vr the Fano factor is found to move from values smaller 
than 1 to values larger than 1 at decreasing values of 
a. We note that for values of a = — 2 < —3/2 in 3D 
(a = —3/2 < —1 in 2D) the sample displays enhanced 
shot noise at increasing voltages, in complete agreement 
with theoretical predictions for d = 2 and d — 3. For 
values of a > —3/2 in 3D (and a > — 1 in 2D) the sys- 
tem displays suppressed shot noise at increasing voltages, 
thus confirming the trend found in previous simulations. □ 
For a — —3/2 in 3D (and a = — 1 in 2D) the system is 
found to display full shot noise within numerical uncer- 
tainty, as expected from theoretical predictions. Finally, 
by comparing the general trend of the Fano factor with 
the I—V characteristics we find that in passing from sup- 
pressed to enhanced shot noise the conductor exhibits a 
substantial nonohmic behavior with a nearly saturation 
current regime at the largest negative a values. 

To detail the evolution from suppressed to enhanced 
shot noise, we have reported in Fig. ^ the value of the 
Fano factor for the 3D and 2D cases at V/Vr = 60 for all 
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values of a considered here. For the purpose of compari- 
son, the figure also reports the theoretical predictions of 
Refs. [ @, for d = 3. The Fano factor is found 
to exhibit an asymmetric behavior with a, showing a 
tendency to saturate, taking suppressed values near to 
l/d, at positive values of a, while increasing systemat- 
ically towards enhanced values well above unity at the 
largest negative values of a. The quantitative agreement 
between MC simulations and existing theories is found 
to be reasonably good for —1.5 < a < 0.5. These the- 
ories are absent or fail completely outside this range of 
values. By contrast, starting from a < —d/2 the simula- 
tions evidence an enhanced shot noise .regime which fully 
validates the theoretical expect at ions .£3 

To investigate the mechanism responsible for the on- 
set of shot-noise enhancement, Fig. H reports the spatial 
profiles of the relevant average quantities (concentration, 
electric field, velocity and scattering time) for an applied 
voltage X/Vt = 80 and for different values of a in a 3D 
system.O From the figure we identify the signature of a 
qualitative change in the different profiles for values of 
a below or about —1. Here, we find the onset of two 
maxima (minima) in the velocity (concentration), more 
pronounced the lower the value of a. This feature de- 
termines the presence of two regions inside the sample 
characterized by two transport regimes: a quasi-ballistic 
one in the region near to the contacts and a diffusive one 
around the center of the structure. Furthermore, a wide 
region of dynamical negative differential mobility shows 
up inside the sample, which in turns moves the structure 
towards a state of electrical instabilitjL,. as well known 
for the analogous case of Gunn diodes.EJ These qualita- 
tive changes are accompanied by a systematic increase of 
the Fano factor, as can be seen in Figs. [I] and |], what 
indicates the onset of a mechanism inducing positive cor- 
relations among the fluctuations. Therefore, the increase 
of the Fano factor stems as a precursor of the fact that 
the system is evolving towards such an electrical instabil- 
ity. The appearance of a negative differential resistance 
in the I — V characteristics (more pronounced in the 2D 
case), see insets of Fig. [j], further supports the present 
interpretation. 

Figure |I] reports the spatial profiles of the same rel- 
evant quantities of Fig. |3| at different applied voltages 
for the 3D case with a = —2. Here the onset of the 
bimodal profiles of carrier concentration and velocity at 
increasing voltages is clearly confirmed, thus completing 
the analysis of the microscopic mechanism responsible of 
the shot-noise enhancement started-with Fig. ||. 

The rjPfWfiflfti MC simulationso'tfB and analytical 
theories0tl'Bli3O have considered applied bias up to val- 
ues high enough to reach the conduction regime known 
as space-jcharge limited conditions and defined by the 
conditionQ L 3> L s ^> Lp or equivalently Vr <C V <C 
V s hot- However, in Sec. [Hl|, starting from suppressed 
shot noise we have predicted the possible achievement 
of full shot noise for bias values satisfying the condition 
V > V s hot, which is in general beyond the space-charge 



limited regime. To test this prediction, we have per- 
formed MC simulations for a 3D system with scattering 
parameters a — 3/2 and r = 8.8 10~ 15 s cV -3 / 2 , at ap- 
plied voltages sufficiently high to satisfy V > V s hot- The 
results of such simulations are reported in Fig. ||, which 
shows the low frequency current spectral density (con- 
tinuous curve and full dots) and the current (multiplied 
by 2q, dotted curve) as a function of bias for the simu- 
lated structure. For the sake of completeness, the inset 
in the same figure shows the corresponding Fano factor. 
The simulations prove that at the highest bias the sample 
indeed displays full shot noise, thus confirming theoret- 
ical expectations. We note that the values of the ap- 
plied voltage needed to achieve full shot-noise conditions 
are much higher than those reported in Fig. |]. Focus- 
ing on the inset, we observe that the crossover between 
thermal and full shot noise does not follow the standard 
monotonic expression 7 = coth(qV/2kBT), but is medi- 
ated by a transition region that displays suppressed shot 
noise down to a minimum value. We remark that for 
the scattering parameters chosen here full shot noise ap- 
pears at high voltages for which the I — V characteristic 
is still far from the current saturation region. There- 
fore, this full shot-noise regime can not be related in any 
way to the presence of current saturation, as it happens 
in other systems like vacuum tubealj or ballistic non- 
degenerate conductors .63611123 Rather, it constitutes an 
intrinsic property of the sample related to the diffusive 
transport regime of carriers (the same phenomenon can 
be observed in inelastic diffusive nondegenerate semicon- 
ductors, see below and Ref. [ ^3|). For the particular 
value of a = 3/2 the minimum value of the Fano fac- 
tor in the anomalous region of the crossover is around 
f/3. For other values of a £ (— f.5,3/2) the value of 
the minimum is found to lay in the range (1,1/3), as 
shown in Fig. g. According to these results, we conclude 
that for sample lengths longer than the contact Debye 
screening length, L > Ld o , and for a > —d/2, a nonde- 
generate elastic diffusive conductor exhibits the follow- 
ing three noise regimes: (i) thermal noise at low bias; 

(ii) suppressed shot noise, with a minimum value of the 
Fano factor determined by the energy dependence of the 
scattering time (a), at an intermediate bias range where 
the sample is under space-charge limited conditions and; 

(iii) full shot noise, independently of a, at the highest 
bias range. 



B. Anomalous crossover between thermal and shot 
noise 

^From the results of the previous section it emerges 
that the shot-noise suppression so found can be more 
generally interpreted in the context of an anomalous 
crossover between thermal and shot noise, as introduced 
in Ref. [ |2^| for the case of macroscopic diffusive con- 
ductors where scattering is inelastic. To shed more light 
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on this subject, we find of interest to carry out a com- 
parison between the noise properties of an elastic non- 
degenerate conductor with those of an inelastic one. To 
this purpose we have performed simulations for a system 
with the same properties of that studied above but in the 
presence of inelastic collisions, using a scattering model 
already developed in Ref . [ [I| . The results are reported 
in Fig. |^, which shows the low frequency current spectral 
density for both elastic and inelastic diffusive conductors 
as a function of the current flowing through the sam- 
ples. For the sake of completeness the inset reports the 
corresponding Fano factor. In both the elastic and in- 
elastic cases, at increasing voltages we can distinguish 
three different noise regimes, corresponding to thermal, 
suppressed and full shot-noise conditions. The degree 
of suppression of shot noise in the intermediate region 
is found to be much more pronounced in the inelastic 
case. This fact is more clearly illustrated in the inset of 
Fig. |^. Remarkably enough, the transition between the 
thermal and the full shot-noise regimes in the inelastic 
case is mediated by a well defined anomalous crossover 
with a cubic dependence of the current spectral density 
upon the current, in agreement with the results of Ref. [ 
. For the elastic case, the crossover displays a power 
dependence with current to an exponent of three halves 
(see Fig. ||), which needs to be further investigated. We 
conclude that the anomalous crossover between thermal 
and full shot noise constitutes a unifying scenario to in- 
terpret shot-noise suppression in nondegenerate diffusive 
conductors. 



V. CONCLUSIONS 



is generated by the simultaneous presence of two distinct 
transport regimes, a ballistic one near to the contacts and 
a diffusive one inside the structure. When the sample 
displays suppressed shot noise the Fano factor exhibits a 
minimum whose value depends on a. This dependence 
is in reasonable agreement with analytical theories de- 
veloped for space-charge limited conditions and d — 3 
in the range —1.5 < a < 0.5. Outside this range ana- 
lytical theories fail completely providing unphysical val- 
ues. Since full shot noise is asymptotically reached at the 
highest bias, shot-noise suppression (and the correspond- 
ing minimum value of the Fano factor) is more generally 
interpreted in terms of an anomalous crossover between 
thermal and full shot noise. We believe that this inves- 
tigation provides a comprehensive understanding of the 
noise properties of nondegenerate elastic diffusive con- 
ductors, and sheds new insight into the nonequilibrium 
noise properties of mesoscopic conductors. 
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We have presented a detailed analysis of the noise prop- 
erties of nondegenerate elastic diffusive conductors. To 
this end, we have performed a complete set of Monte 
Carlo simulations on the basis of the indications pro- 
vided by theoretical predictions. The simulations have 
confirmed the existence of enhanced shot noise when the 
exponent a of the energy dependence of the relaxation 
time satisfies the condition a < —d/2, with d the di- 
mensionality of the system in momentum space. Fur- 
thermore, the presence of full shot noise is confirmed un- 
der different conditions, namely: (i) when the length of 
the sample is shorter than the contact Debye length, (ii) 
when a > —d/2 and the applied bias is high enough so 
that the transit time becomes the shortest time scale of 
the system, and (iii) when a = —d/2. The presence of 
enhanced shot noise for a < —d/2 is found to coincide 
with the analytical condition that the differential conduc- 
tivity in energy space becomes negative. The presence of 
super-Poissonian values of the Fano factor is here inter- 
preted as precursor of an electrical instability driven by a 
region of dynamical negative differential mobility inside 
the sample. This instability is associated with the bi- 
modal spatial profile of the carrier average velocity, which 
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